Astronomy & Astrophysics manuscript no. 44Tau2009 


© ESO 2009 


November 23, 2009 





A 6 Scuti star in the post-MS contraction phase: 44 Tauri 

P. Lenz 1 , A. A. Pamyatnykh 1,2 ' 3 , T. Zdravkov 2 , and M. Breger 1 



o 
o 

> 

o 
m 

(N 
P? 

43 

Or 

6 

o3 



1 Institut fur Astronomie, University of Vienna, Turkenschanzstrasse 17, A- 11 80 Vienna, Austria 
e-mail: patrick . lenz@univie . ac . at 

2 Copernicus Astronomical Centre, Bartycka 18, 00-716 Warsaw, Poland 

3 Institute of Astronomy, Russian Academy of Sciences, Pyatnitskaya Str 48, 109017 Moscow, Russia 



> 
ON 
CO 
CO 

^1- 



o 



Received date; accepted date 



ABSTRACT 



Context. The evolutionary stage of the 6 Scuti star 44 Tau has been unclear. Recent asteroseismic studies have claimed models on the 
main sequence, as well as in the expansion phase of the post-main sequence evolution. However, these models could not reproduce all 
of the observed frequencies, the mode instability range, and the fundamental stellar parameters simultaneously. A recent photometric 
study has increased the number of detected independent modes in 44 Tau to 15, and a newly found gravity mode at 5.30 cd~' extends 
the observed frequency range. 

Aims. One of the possible evolutionary stages of 44 Tau has not yet been considered: the overall contraction phase after the main 
sequence. We computed asteroseismic models to examine whether models in this evolutionary stage provide a better fit of the observed 
frequency spectrum. 

Methods. We used Dziembowski's pulsation code to compute nonadiabatic frequencies of radial and nonradial modes. Observation 
of two radial modes and an avoided crossing of dipole modes put strong constraints on the models. A two-parametric overshooting 
routine is utilized to determine the efficiency of element mixing in the overshoot layer above the convective core. 
Results. We find that pulsation models in the post-MS contraction phase successfully reproduce the observed frequency range, as 
well as the frequency values of all individual radial and nonradial modes. The theoretical frequencies of the mixed modes at 7.79 cd~' 
and 9.58 cd~' are in better agreement with the observations if efficient element mixing in a small overshoot layer is assumed. 

Key words, stars: variables: 6 Set - stars: oscillations - stars: individual: 44 Tau 



1. Introduction 

The 6 Scuti stars are a class of pulsators exhibiting low-order 
acoustic and gravity modes driven by the opacity mechanism 
acting in the Hell ionization zone. In evolved S Scuti stars mixed 
modes, i.e., modes with acoustic behaviour in the envelope and 
gravity behaviour in the interior, are common. Since such modes 
probe different layers inside a star than purely acoustic modes, 
the 5 Scuti stars are excellent targets for asteroseismology. 

44 Tau is a 6 Scuti pulsator of spectral type F2 IV. Its ex- 
ceptionally slow intrinsic rotation of V rot = 3+2 km s (Zima 
et al. 2007) is much lower than the average value for 6 Scuti 
stars. For this group of pulsators, projected rotational velocities, 
vsim, around 100 km s _1 and higher are not unusual. This makes 
44 Tau an interesting and relatively simple target for asteroseis- 
mology, because the effects of rotation are weak and do not com- 
plicate mode identification and stellar modelling. 

In recent years several ground-based observing campaigns 
of the Delta Scuti Network (DSN) were organized to retrieve 
photometric and spectroscopic data of 44 Tau. The first com- 
prehensive frequency analysis of its photometric light variations 
was published by Antoci et al. (120071 1. This paper also summa- 
rizes all previous observations of the star. With two additional 
seasons of extensive Stromgren vy photometry, Breger & Lenz 
(2008) could increase the number of detected frequencies to 49, 
of which 15 are independent. 

Garrido et al. (2007) have determined the spherical degree, £, 
for the 10 dominant modes by means of photometric amplitude 
ratios and phase differences. Their results were later confirmed 
by Lenz et al. (2008). In 2004 a spectroscopic campaign was or- 



ganized simultaneously to the photometric observing run. The 
results of the line-profile variation analyses and the determina- 
tion of photospheric element abundances are discussed in detail 
by Zima et al. (2007). For 8 modes, the azimuthal order, m, could 
be determined. Moreover, no indications of a global magnetic 
field and no anomalous element abundances were found. 

The spectroscopically determined logg value of 3.6 + 0.1 
does not allow for unambiguous determination of the evolution- 
ary status of 44 Tau. In recent studies, both main sequence (here- 
after MS) and post-MS models in the expansion phase have been 
considered (Garrido et al. 2007, Lenz et al. 2008). The results by 
Lenz et al. (2008) showed that these models fail to reproduce all 
observable parameters satisfactorily (individual frequencies, ob- 
served frequency range, fundamental stellar parameters). 

The effective temperature and luminosity of MS models is 
significantly lower than the values derived from photometry and 
Hipparcos data. Moreover, the positions of predicted 1=2 modes 
do not match the observed frequencies (see Fig. 8 in Lenz et al. 
2008). 

Evolved models in the hydrogen-shell burning stage repro- 
duce the fundamental stellar parameters well, but theory pre- 
dicts many more frequencies than actually observed. While par- 
tial mode trapping may explain the position of the €=\ modes, 
this mechanism should not be very effective for the observed 
1=2 modes. Moreover, the frequency range of unstable modes 
is shifted to higher frequencies with respect to MS models. 
Therefore, the mode at 6.34 cd 1 and the newly found mode at 
5.30 cd 1 are predicted to be stable by models in the post-MS 
expansion phase. 
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Table 1. Observed frequencies and mode identifications. 
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We reexamined asteroseismic models of 44 Tau based on the 
latest frequency solution of Breger & Lenz (2008) and inves- 
tigated a hitherto neglected possibility: a model in the overall 
contraction phase after the main sequence. 

During the main sequence evolution, a 6 Scuti star trans- 
forms hydrogen into helium in its convective core. The amount 
of hydrogen in the core gets very low towards the end of this 
phase. The star increases its central temperature by an overall 
contraction to keep the energy production from hydrogen burn- 
ing efficient. At the end of this overall contraction stage, hydro- 
gen in the core is fully depleted and the convective core dis- 
appearfl A nuclear hydrogen-burning shell is established out- 
side the hydrogen-exhausted core. The stellar envelope expands 
again, which leads to a decrease of the effective temperature. 

Finding a 5 Scuti star in the contraction phase after the 
TAMS is not unlikely. If we consider the same range of effective 
temperatures, the evolution during the overall contraction phase 
is approximately 10 times faster compared to the MS evolution. 
In the post-MS expansion phase, the evolution is approximately 
16 times faster than on the MS. The asteroseismic inferences 
from post-MS contraction models are presented in this paper. 

The paper is organized as follows. In Sect. 2 we briefly sum- 
marize the observed parameters of 44 Tau. In Sect. 3 we discuss 
the computation of our seismic models. The diagnostic power 
of mixed modes is examined in Sect. 4. Finally, in Sect. 5 we 
present predictions of post-MS contraction models. 

2. Observed frequencies and mode identification 

All observed independent frequencies of 44 Tau are given in 
Table [TJ along with their mode identifications, if available. The 
given A v amplitudes are mean Stromgren y amplitudes derived 
from the five observing seasons between 2000 and 2006. The 
amplitudes for the first moments are taken from Zima et al. 
(2007). The spherical degrees of the modes were determined by 
Garrido et al. (2007) and Lenz et al. (2008). We reinvestigated 



1 Of course, convective overshooting also disappears, but we keep 
the term 'model with overshooting' or 'the overshooting case' for later 
evolutionary stages to distinguish this case from models constructed 
without convective overshooting. 



these mode identifications with the new full frequency solution 
of Breger & Lenz (2008) based on photometric amplitude ra- 
tios and phase differences and using the method by Daszynska- 
Daszkiewicz et al. (2003 2005 ). The previous results were con- 
firmed. The new mean amplitude ratios and phase differences 
between the Stromgren v and y bands are listed in TableQ] These 
average values were computed from the annual solutions of the 
five observing seasons and weighted by the annual signal-to- 
noise ratio of each mode. Taking the observed amplitude modu- 
lation in 44 Tau into account, this weighting scheme is more ad- 
equate than the one used in Lenz et al. (2008). For some modes, 
in particular for fi), the uncertainties in A V /A Y and <p v - <f> y could 
be reduced significantly. 

The spherical degree can be uniquely determined for the 10 
dominant modes. The observed mean phase difference of fu 
puts this mode midway between the theoretical values of I — 1 
and 2 modes (see Fig. 5 in Lenz et al. 120081 ). which prevents 
an unambiguous identification. The uncertainties in the phase 
differences of the low-amplitude modes fn-fis are too large to 
derive reliable results. 

The mode positions in the diagnostic diagrams for mode 
identification depend on the efficiency of convection. By fit- 
ting the predicted mode positions to the observed values, Lenz 
et al. (2008) determined the mixing-length parameter, q^mlt, to 
be < 0.2 in 44 Tau. This is also confirmed by the new data. 
Consequently, convection in the envelope of 44 Tau is not very 
efficient. 

The azimuthal order of eight modes was derived by Zima et 
al. d2007 ) using the moment method. Most modes in 44 Tau are 
axisymmetric with exception of the two prograde modes fi and 
fi. No rotational splitting has been observed. The close pair f(, 
and ft) at 9.56 and 9.58 cd~' does not represent a part of a rota- 
tionally split multiplet, because the determined spherical degrees 
differ. The inclination angle of 44 Tau could be constrained to 60 
+ 25°. 

The effective temperature of 44 Tau is r eff = 6900 + 100 K, 
logL = 1.305 ± 0.065, and logg = 3.6 ± 0.1. Due to the use of 
a more recent estimate for the bolometric correction the given 
mean luminosity is lower than the value in Lenz et al. (2008). 
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Table 2. Parameters of post-MS contraction models with a good fit of the 15 frequencies observed in 44 Tau. 
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3. Seismic models in the post-MS contraction 
phase 

In this study we used the Warsaw-New Jersey evolutionary 
code, which is a descendant of the Paczyriski (1970) code. In 
this code, envelope convection is described with the standard 
mixing-length theory of convection. The evolutionary models 
are computed assuming uniform rotation and conservation of 
global angular momentum during evolution from the ZAMS. 

Aside from the traditional description of overshooting from 
the convective core, we also applied a new two-parametric de- 
scription that allows for partial element mixing in the overshoot- 
ing region (Dziembowski & Pamyatnykh 2008). It allows us 
to consider different profiles of the hydrogen abundance inside 
the partly mixed region just above the convective core. One of 
the parameters, a ov , defines the overshooting distance above the 
convective core. The second parameter, w, adjusts the profile of 
the hydrogen abundance between the convective core and the up- 
per boundary of the overshooting region. It therefore character- 
izes the efficiency of mixing in this region. The hydrogen mass 
fraction in a stellar layer, X, is defined as a function of fractional 
mass, q: 

X = X c + (q-q c y[a + b(q-q c )]. (1) 

The coefficients a and b are determined from w and a ov . The 
subscript c refers to the central values of the given quantities. An 
infinitely high value of w corresponds to the traditional standard 
treatment for overshooting. 

To compute the nonadiabatic radial and nonradial frequen- 
cies, we used the most recent version of Dziembowski's pulsa- 
tion code d 19771 1. This code relies on the frozen convective flux 
approximation. The effects of rotation on the frequencies were 
treated up to second order following Dziembowski & Goode 
(fT992l . 

In this study we used two sources of the Rosseland mean 
opacities: the OPAL opacity tables (Iglesias & Rogers 1996) 
and the OP opacities (Seaton 2005). At low temperatures, these 
tables are supplemented by the Ferguson & Alexander (2005 ) 
molecular and grain opacities. 

The abundances of some photospheric elements in 44 Tau 
were determined by Zima et al. (2007). Within the given uncer- 
tainties, the abundances of 44 Tau are close to the solar values 
by Grevesse & Noels 019931 hereafter GN93) and Asplund et al. 
( 120041 120531 hereafter A04). Since only abundances of 1 1 ele- 
ments were measured, we did not utilize a special element mix- 



ture for 44 Tau, but examined pulsation models with both solar 
element mixtures. 

3.1. Results with OPAL opacities 

Two radial modes are observed in 44 Tau. As discussed in Lenz 
et al. (2008), this allows us to utilize Petersen diagrams (Petersen 
& J0rgensen 1972). In these diagrams, the period ratio between 
two radial modes is plotted against the period of the fundamen- 
tal mode. The period ratio is sensitive not only to the adopted 
metallicity but also to rotation as stated by Suarez et al. (f2006 ). 
Since these quantities are determined well by observations, we 
can use Petersen diagrams to find the optimum mass of a model. 

Another important observation puts tight constraints on the 
models: the small frequency separation of the two £ = 1 modes 
at 9.11 and 9.56 cd~' indicates that these modes undergo an 
avoided crossing (Aizenman et al. 1977). Since both modes are 
of mixed character, they are also sensitive to the extent of over- 
shooting from the convective core. 

To find the optimum fit for the dipole modes, a series of 
computations with different overshooting distances, a ov , and dif- 
ferent efficiencies of partial element mixing in the overshoot 
layer (w=2.0 and 8.0) were made. Models that exhibited a good 
fit of all dipole modes also showed very good agreement for 
the quadrupole modes. For models obtained with X=0.70, the 
frequency separation between the dipole modes at 9.11 and 
9.56 cd -1 is predicted to be slightly larger than observed. An 
increase in the initial hydrogen mass fraction, X, helps to de- 
crease the predicted frequency separation to the observed values 
and provides the best agreement for all observed 15 frequencies. 

The best pulsation models for OPAL opacities and the GN93 
and A04 element mixture are given in Table [2] All models were 
computed with a mixing length parameter <?mlt=0.2. Model 4 
was computed with the traditional description of overshooting. 
Starting with an initial rotation rate of 3.5 km s _1 , the tabulated 
evolved models exhibit rotation rates close to 3 km s _1 . The ef- 
fects of rotation were also taken into account in the calculation 
of the pulsation frequencies. 

A comparison between predicted and observed frequencies 
for Model 4 in Table [2] is shown in Fig. Q] It can be seen that 
this pulsation model allows all observed modes to be explained 
solely with ( < 2 modes. In some cases modes with {-3 or 4 
are also close to observed values, but photometric cancellation 
effects clearly favour low degree modes. 
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Fig. 1. Comparison between observed frequencies (vertical 
lines) and predicted unstable frequencies (circles) for a post-MS 
contraction model (Model 4 in Table 13. If the spherical degree 
of a mode was not observationally determined, a full vertical line 
is shown. Predicted frequencies that match observed frequencies 
are marked as filled circles. The rotational splitting is smaller 
than the size of the symbols. 



To describe the goodness of the fit, we used a dimensionless 
merit function similar to the one used in Brassard et al. (2001): 
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where v^ bs and v| nodel are the observed and corresponding theo- 
retical frequency, respectively. The weighting factor, cr 1 , is de- 
fined as the width of the observed frequency window of 5- 
13 cd~' divided by the number of theoretically predicted modes 
within this frequency window. Contrary to Brassard et al. (2001) 
we determine a weighting factor for each spherical degree sepa- 
rately to take the different mode densities into account. For the 
models shown in Tabled cr' for t — modes amounts to 8/3, 
for ( — 1 modes 8/7 and for t — 2 modes 8/11 cd~'. The re- 
sults are shown in the last column of Table [2] The lowest x 1 
value provides the best agreement between the model and the 
observations. If the azimuthal order of a mode was not uniquely 
determined, the theoretical m = mode was used to compute^ 2 . 
If the m-component closest to the observed frequency is used in- 
stead, the x 1 values are slightly lower but exhibit the same gen- 
eral trends. 

The position of OPAL models in the HR diagram is shown in 
Fig.|2]for models obtained with the GN93 and the A04 element 
mixture. The effective temperatures of these pulsation models 
are somewhat cooler than the values derived from photometric 
and spectroscopic measurements. 

The pulsation models obtained with the A04 mixture are 
closer to the TAMS. A 15-frequency fit with w=2.0 is not pos- 
sible, because the required overshooting distance would be less 
than possible for models in the overall contraction phase. The 
main reason for the differences between pulsation models ob- 
tained with the GN93 and the A04 mixture is the slightly differ- 
ent opacity in the overshoot layer that affects the size of the con- 
vective core. The uncertainties in the element abundances (and 
opacities) in 44 Tau therefore lead to an uncertainty in the deter- 
mination of the overshooting distance, a ov . 
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Fig. 3. Instability parameter, rj, for Model 3. Positive values in- 
dicate mode driving, modes with negative rj are damped. 



The predicted frequency range of unstable modes agrees well 
with the observed frequency range as shown in Fig.[3]for one of 
the models that can be considered representative of all OPAL 
models. Additional unstable modes are predicted between 13 
and 17 cd . The mode instability at these high frequencies is 
sensitive to the efficiency of convection. If ctmlt is decreased 
from 0.2 to lower values, the highest frequencies become stable. 
In our study we relied on the standard mixing-length theory of 
convection and the frozen flux approximation. Time-dependent 
convection models, such as those used in some other recent stud- 
ies about 6 Scuti stars (e.g., Montalban & Dupret 2007, Dupret 
et al. 120051 1. would yield a more accurate determination of the 
high-frequency border of instability. 

3.2. Results with OP opacities 



Lenz et al. (120071) found that post-MS expansion models of 
44 Tau constructed with OP opacities have significantly lower 
temperature and luminosity values than observed. This problem 
persists for post-MS contraction models. Montalban & Miglio 
(2008) explain this by differences of 10% between OPAL and 
OP opacities at temperatures around log T = 6.05. This corre- 
sponds to the temperature region in which the period ratio is a 
sensitive probe as can be seen in the two uppermost panels in 
Fig.0 

The diagram shows the kinetic energy density of the radial 
fundamental and first overtone mode inside the star. Pulsation 
modes are sensitive to the conditions in temperature regions in 
which the kinetic energy density is high. Because of the position 
of its node, the radial first overtone mode probes the temperature 
region between log T = 4.5 and 7.0 with different weights than 
the radial fundamental mode. 

The differences between OPAL and OP opacities are shown 
in the middle panel of Fig. |4] The lower panel shows the location 
of opacity bumps inside the star. We artificially enhanced the OP 
opacities by up to 15% around log T = 6.0 and tested the impact 
on models in the contraction phase after the main sequence. 

The results are shown in Fig. [5] The corresponding model pa- 
rameters can be found in Table|2] Contrary to models constructed 
with standard OP opacities, models obtained with enhanced OP 
opacities are closer to the photometric error box in the HR dia- 
gram. This may indicate that the Rosseland mean opacities of the 
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(a) Models constructed with the GN93 element mixture. 
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(b) Models constructed with the A04 element mixture. 
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Fig. 2. HR diagram with evolutionary tracks for models with a good fit of all 15 observed modes. The models were constructed with 
OPAL opacities and two different solar element mixtures: the GN93 element mixture (a) and the A04 element mixture (b). 




Fig. 4. Radial fundamental and first overtone mode as probes 
for opacities. Upper panel: kinetic energy density inside the 
star. Middle panel: comparison between OPAL and OP opacities 
for the same model. Lower panel: absolute values of Rosseland 
mean opacities from OPAL and OP data. 



OP tables are currently underestimated in the temperature region 
around log T = 6.0. 

For models constructed with standard OP opacities, a fit of 
all 15 frequencies is only possible with the GN93 element mix- 
ture. Adopting OP instead of OPAL opacities shifts the mod- 
els that fit all 15 frequencies closer to the TAMS (e.g., com- 
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Fig. 5. HR diagram with evolutionary tracks of two models con- 
structed with enhanced and standard OP opacities. Both models 
provide a good 15 -frequency fit. 

pare the models in Fig. [2^ with Model 7 in Fig. [5j. The models 
constructed with OPAL A04 are already situated at the TAMS 
(Fig. I2J5). To fit the avoided crossing of I = 1 modes and the 
I — 2 modes in standard OP A04 models, a m would have to be 
reduced more than possible for a model in the overall contraction 
phase. However, the use of enhanced OP opacities in combina- 
tion with the A04 element mixture allows for a 15-frequency fit. 
The corresponding models are listed in Table [2] 

4. Mixed modes as probes of the stellar core 

The partial mixing processes at the convective core boundary are 
currently not fully understood. During the hydrogen core burn- 
ing phase, we expect a region around the convective core that 
is partially mixed due to several mechanisms: (i) overshooting 
from the convective core, (ii) rotationally induced mixing (the 
core rotates faster than the envelope leading to additional mix- 
ing). 
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Unfortunately, until now no rotational splittings could be de- 
tected in 44 Tau. Therefore, we cannot derive the profile of dif- 
ferential rotation that would be necessary to disentangle the ef- 
fects of rotation in element mixing. 

An excellent probe of the size of the overshoot layer is the g\ 
mode for I > as stated by Dziembowski & Pamyatnykh (1991). 
This mode is partially trapped in the overshoot region, so its fre- 
quency has high diagnostic value. However, in 44 Tau models in 
the contraction phase after the TAMS, this mode already moved 
to higher frequencies outside the observed frequency range. 

However, we observe many mixed modes in 44 Tau which 
are also sensitive to the conditions in the chemically inhomoge- 
neous overshoot layer above the convective core. Table [3] lists 
the fraction of oscillation kinetic energy confined in the g-mode 
cavity, £kin,g/£km, for all predicted unstable modes in Model 3. 
Since the given values of Zikin.g/^km are very similar for the dif- 
ferent OPAL and OP models, they can be considered represen- 
tative of all post-MS contraction models. Moreover, the differ- 
ences between calculated and observed frequencies are given. If 
the azimuthal order of an observed mode is unknown, the fre- 
quency difference to the theoretical m=0 modes is given and the 
value is enclosed by brackets. The uncertainty in frequency due 
to unknown azimuthal order may be as high as 0.05 cd for 
I — 2 modes. 
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Fig. 6. Comparison of the theoretical frequencies of two pulsa- 
tion models with different hydrogen profile modelling, w=8.0 
(filled circles) and w=2.0 (open circles). Modes that are most 
affected are marked with a dotted circle. The parameters of the 
models are given in Table[2](see Models 2 and 3). 



Table 3. Theoretical frequencies and mode properties of unsta- 
ble modes predicted by Model 3. 



Mode 


"model 


^observed 


^observed ^model 


^kin,g / £kin 


t 


ID 


[cd- 1 ] 


[cd- 1 ] 


[cd- 1 ] 







F 


6.8980 


6.8980 


0.0000 


0.0 





1H 


8.9607 


8.9606 


-0.0001 


0.0 





2H 


11.20 






0.0 





3H 


13.48 






0.0 




g6 


5.44 






0.96 




g5 


6.29 






0.92 




Pi 


7.0342 


7.0060 


-0.0282 


0.13 




g4 


7.7915 


7.7897 


(-0.0018) 


0.87 




P2 


9.1099 


9.1174 


0.0075 


0.35 




g3 


9.5748 


9.5611 


(-0.0137) 


0.61 




P3 


11.5526 


11.5196 


-0.0330 


0.04 




P4 


13.74 






0.19 




g2 


14.51 






0.75 


2 


gin 


5.3375 


5.3047 


(-0.0328) 


0.96 


2 


g9 


5.78 






0.93 


2 


g8 


6.3076 


6.3390 


(0.0314) 


0.83 


2 


g7 


6.7795 


6.7955 


0.0160 


0.68 


2 


f 


7.3023 


7.3031 


0.0008 


0.67 


2 


g(> 


7.98 






0.74 


2 


Pi 


8.6401 


8.6391 


-0.0010 


0.64 


2 


g5 


9.5584 


9.5828 


(0.0244) 


0.62 


2 


P2 


10.47 






0.57 


2 


g4 


11.3129 


11.2947 


(-0.0182) 


0.60 


2 


P3 


12.7043 


12.6915 


(-0.0128) 


0.34 


2 


g3 


13.92 






0.04 


2 


P4 


15.17 






0.35 



of other dipole modes in 44 Tau. Consequently, this and other 
mixed modes are sensitive to the shape of the hydrogen profile 
in the partially mixed layer above the convective core. 

We computed pulsation models using two different assump- 
tions for element mixing in the overshoot layer: weak element 
mixing (w=2.0) and efficient element mixing (w=8.0). For com- 
parisons we also calculated models with the standard description 
of overshooting, which assumes no element mixing in the over- 
shooting region. 

In Fig. [6] we compare the predicted frequency spectra for 
Models 2 and 3 (obtained with w=2.0 and w=8.0, respectively). 
The corresponding overshooting distances were adjusted to val- 
ues that provide the best fit of the observed frequencies. As can 
be seen, the frequencies of some modes are sensitive to the dif- 
ferences in the overshooting region. The frequency of the g4 
mode at 7.79 cd -1 shifts most significantly. The frequency dif- 
ference Av=0.1 cd 1 is larger than the expected maximum fre- 
quency separation between components of a rotational splitting. 
Consequently, our results suggest preference for more efficient 
mixing in a small partially mixed region (w=8.0) over less effi- 
cient mixing in a larger partially mixed region (w=2.0). 

The corresponding profile of the hydrogen abundance for 
the discussed pulsation models is given in the upper panel of 
Fig. [7] The different hydrogen profiles lead to a modification 
of the Brunt-Vaisala' frequency, as shown in the panel in the 
middle of the diagram. In the lower panel the eigenfunction of 
the relative Lagrangian pressure variation is given for the mixed 
I = 1 modes: gn at 7.79 cd 1 and p\ at 7.01 cd -1 . Due to its 
larger amplitude in the partially mixed region the gravity mode 
is very sensitive to the conditions in this region. The change in 
the Brunt-Vaisala frequency due to the different chemical profile 
slightly adjusts the size of the g-mode cavity and, therefore, also 
influences the values of the frequencies. 



The mixed C-l mode at 7.79 cd has strong g-mode char- 
acteristics. Almost 87% of its kinetic energy is confined in the g- 
mode cavity in the stellar interior. This also explains why the ob- 
served amplitude of this mode is significantly smaller than that 



5. Predictions of the post-MS contraction model 

Post-MS contraction models of 44 Tau predict a few additional 
dipole and quadrupole modes that may be found by future ob- 
servations. After fitting the observed 15 modes, the frequencies 
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Fig. 7. Impact of different efficiency of element mixing in the 
overshoot layer on the profile of the hydrogen abundance (up- 
per panel) and on the Brunt-Vaisala frequency (middle panel). 
Lower panel: relative Lagrangian pressure perturbation for two 
mixed £=l modes, the mainly acoustic mode at 7.01 cd l and 
the mainly gravity-type mode at 7.79 cd~' . 



of the remaining modes no longer strongly depend on the input 
parameters. However, small deviations may be found, especially 
for g-modes that will provide additional constraints on mixing 
in the overshoot layer. The expected frequencies of hitherto not 
observed modes can be found in Table [3] The detection of these 
modes would help refining our models. 

A reexamination of our photometric data indeed shows a 
prominent peak at the predicted position of the second radial 
overtone (1 1.198 cd" 1 ). With an SNR of 2.72, this peak is below 
the significance limit of 4.0 but may be confirmed with new addi- 
tional data. A reexamination of radial velocity data did not show 
any reliable prominent peaks at the predicted position of modes. 
Many of the predicted modes that lack an observational coun- 
terpart are gravity modes for which we only expect low ampli- 
tudes. Therefore, more accurate data are needed to extract these 
frequencies. 

Moreover, the model predicts the spherical degree of some of 
the observed modes: € = 1 for 7.79 cd" 1 and i = 2 for 5.30, 8.64, 
11.30, 12.69 cd~'. These predictions also need to be confirmed 
by high-resolution spectroscopic data. 

The computed evolutionary relative period changes, 
(l/P)dP/dt, are approximately 1 • 10~ 8 yr _1 . Such changes are 
too small to be measured with the 5 year time-base of our data. 



6. Conclusions 

We constructed seismic models of the 5 Scuti star 44 Tau in the 
overall contraction phase after the main sequence. The two ob- 
served radial modes and an avoided crossing of a pair of dipole 
modes put strong constraints on pulsation models. Unlike mod- 
els on the main sequence or in the post-MS expansion phase, 
these models successfully reproduce all 15 independent frequen- 
cies observed in 44 Tau and predict all detected modes unstable. 
This makes 44 Tau an example for determining the evolutionary 
stage by means of asteroseismology. 

The observed frequency spectrum (obtained from ground- 
based data) can be solely explained with modes with spherical 
degrees of 0, 1 and 2. This result is what we expect for the given 
accuracy of the data. 

The frequency spectrum of 44 Tau consists mainly of modes 
of mixed character, i.e., modes with acoustic behaviour in the 
envelope and gravity behaviour in the interior. These modes al- 
low examination of the efficiency of partial element mixing in 
the overshoot layer around the stellar core. We find that the the- 
oretical frequencies of several modes (e.g. the dipole g4 mode at 
7.79 cd~' and the quadrupole mode gs at 9.58 cd~') are in better 
agreement with efficient mixing in a thin overshoot layer than 
with less efficient mixing in an extended overshoot layer. 

The main source of uncertainty in the asteroseismic mod- 
elling of 44 Tau are the stellar opacities. Determination of the 
size of the overshooting region above the convective core, there- 
fore, depends on the source of opacity data and, to a smaller 
extent, on the choice of the element mixture. Pulsation models 
in the post-MS contraction phase constructed with standard OP 
data are more than two standard deviations outside the photo- 
metric error box in the HR diagram. Following a suggestion by 
Montalban & Miglio (2008), we computed a model with arti- 
ficially enhanced OP opacities around the temperature region 
log T = 6.0. Such a model indeed agrees better with the fun- 
damental parameters of 44 Tau derived from photometry and 
Hipparcos data. 

The pulsation models presented in this paper not only pro- 
vide an excellent fit of all observed modes but also allow sev- 
eral predictions, which may be confirmed with additional ob- 
servational data. A few of the theoretically predicted modes are 
not observed. The question arises whether the amplitudes of 
these modes are below the noise threshold of the current data 
or whether they are even excited at all. Moreover, we predict the 
spherical degree of already observed modes for which our mode 
identification techniques did not succeed because the uncertain- 
ties in amplitudes and phases are too large. The excellent fit of 
the individual frequencies is very promising, and we are con- 
fident that additional data will lead to even more asteroseismic 
inferences about this star. 
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